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Abstract 

We rigorously investigate the quantum dissipative dynamics of a ratchet system de- 
scribed by a periodic potential model based on the Caldeira-Leggett Hamiltonian with a 
biharmonic force. In this model, we use the reduced hierarchy equations of motion in the 
Wigner space representation. These equations represent a generalization of the Gaussian- 
Markovian quantum Fokker-Planck equation introduced by Tanimura and Wolynes (1991), 
which was formulated to study non-Markovian and non-perturbative thermal effects at 
finite temperature. This formalism allows us to treat both the classical limit and the tun- 
neling regimes, and it is helpful for identifying purely quantum mechanical effects through 
the time evolution of the Wigner distribution. We carried out extensive calculations of the 
classical and quantum currents for various temperatures, coupling strengths, and barrier 
heights. Our results reveal that at low temperature, while the quantum current is larger 
than the classical current in the case of a high barrier, the opposite is true in the case of a 
low barrier. We find that this behavior results from the fact that the tunneling enhances 
the current in the case of a high barrier, while it suppresses the current in the case of a 
low barrier. This is because the effect of the ratchet potential is weak in the case of a 
low barrier, due to the large dispersion of the distribution introduced by tunneling. This 
causes the spatio-temporal asymmetry, which is necessary for ratchet current, to be weak, 
and as a result, the net current is suppressed. 

1 Introduction 

A system that is able to rectify thermal or mechanical fluctuations through a periodic potential 
is called a ratchet. P,l2l[3llll|5l|6l[71|8] With recent advances in microscopy and nano-technology, 
several new mechanisms of rectification of fluctuations resulting in transport have been identi- 
fied through theoretical works P [IHl [IIl[I21[I31[Il[ISl[ISl[ni[ISl[IS]and experimental works [201 
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iniESlESlEllESlEniETlEHlEHlEO^ biology, physics, and chemistry. Well-known examples in- 
clude asymmetric quantum dots [201 EH [22] , vortices in superconductors |23 [ [2H 125] . cold atoms 
in asymmetric optical lattices [221 123 12H] , and molecular motors in biological systems. [221 EOl EI] 
A ratchet system can produce a directed current from fluctuations only when the system is in 
a non-equilibrium state, because the extraction of work from unbiased fluctuations is not al- 
lowed by the second law of thermodynamics. Moreover, in order to realize the ratchet effects, 
it is necessary to break the spatio-temporal symmetries, because otherwise the contributions 
from positive and negative currents would cancel. [321 |33l EH ES, EE] Quantum mechanical 
effects may also play a role, because ratchet systems must be of microscopic size to rectify the 
fluctuating motion. [371EHlEnilinillIlll21ll31ll3llSllinillZ] 

A Fokker-Planck approach [18] and a Langevin approach [19] have been used in classical 
studies of ratchets, while varieties of equation of motion approaches have been employed for 
quantum mechanical studies. [371EHlE9lllQll3llll21ll3llllll5lll6lll7] In quantum mechanical 
case, dissipative systems are commonly modeled as potential systems coupled to heat-bath 
degrees of freedom at finite temperature. This coupling gives rise to thermal fluctuations 
and dissipation that drive the systems toward the thermal equilibrium state. The heat-bath 
degrees of freedom are then reduced using such methods as the projection operator method 
and the path integral method, for example. The quantum Langevin equation [37] and the 
quantum Fokker-Planck equation [4:2j have been used for the purpose of understanding the 
quantum aspects of ratchet dynamics in a heat bath. Although these equations are analogous 
to the classical kinetic equations, which have proved to be useful in the treatment of classical 
transport problems, such equations cannot be derived in a quantum mechanical framework 
without significant approximations and/or assumptions. For example, in dealing with the 
quantum Langevin equation expressed in operator form, it is generally assumed that the anti- 
symmetric correlation function of the noise is positive, but this is valid only for slow non- 
Markovian modulation at high temperature, as we demonstrate in Sec. 2. Similarly, the 
quantum Fokker-Planck equation can be derived from the Caldeira-Leggett Hamiltonian under 
a Markovian approximation, [^Tl IS^ 15 ^ IM f 155] but in order for this to be possible, the heat bath 
must be at a sufficiently high temperature, in which case quantum tunneling processes play a 
minor role. The quantum master equation expressed in terms of Floquet states provides a more 
rigorous treatment of quantum dissipative dynamics than the methods mentioned above. [lU] 
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but it can only be applied to systems possessing weak interactions. [50] The methods discussed 
to this point represent the main approaches used in the treatment of ratchet systems carried out 
to this time. However, as we have pointed out, none of them provide fully quantum mechanical 
descriptions of broad validity. 

In this paper, we employ the reduced hierarchy equations of motion (HEOM) in the Wigner 
space representation, [56] which are a generalization of the Gaussian-Markovian quantum Fokker- 
Planck equation introduced by Tanimura and Wolynes.[Fn \5E\ to study non-Markovian and 
non-perturbative thermal effects at finite temperature. [591 EDI EH ESI |63] Because the HEOM 
are derived from the full system-bath Hamiltonian with no approximation, the entire sys- 
tem approaches a thermal equilibrium state at finite temperature when no external perturba- 
tion is applied. [561 El] This means that the second law of thermodynamics holds within the 
HEOM approach and hence there can be no finite current in the absence of a driving force, 
which is the essential behavior to investigate the ratchet rectification of thermal fluctuations. 
The HEOM have been used to study a variety of phenomena and systems, including exciton 
transfer [651 EEl EH ESI EH], electron transfer [TOl [TH CS], quantum information [711 175] . and 
resonant tunneling diodes. [621 [63] The HEOM are ideal for studying quantum transport sys- 
tems when implemented using the Wigner representation, because they allow the treatment of 
continuous systems utilizing open boundary conditions and periodic boundary conditions [76]. 
Elucidation of the dynamical behavior of the system through the time evolution of the Wigner 
distribution functions is also possible. The classical hierarchy equations of motion can be 
obtained easily by taking the classical limit of the HEOM and, utilizing the Wigner representa- 
tion, we can easily compare quantum and classical distribution functions. [5S1 EI] In this paper, 
considering both classical and quantum mechanical cases, we report the results of numerical cal- 
culations of the ratchet current, which is the induced net current resulting from ratchet effects, 
in a dissipative environment for various temperatures and system-bath coupling strengths. We 
then clarify the roles of fluctuations, dissipation and the biharmonic force in both the classical 
and quantum cases. 

The organization of the paper is as follows. In Sec. 2, we introduce the HEOM approach 
in the Wigner representation. In Sec. 3, we discuss ratchet systems from the point of view 
of space-time symmetries. In Sec. 4, the numerical results obtained for the ratchet current in 
classical and quantum mechanical cases are presented. Section. 5 is devoted to conclusions. 
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2 Reduced Hierarchy Equations of Motion Formalism 



Ratchet systems are often modeled by a Brownian particle in a periodic potential under an ac 
or dc driving force. We take this approach and employ a model based on the Caldeira-Leggett 
(or Brownian) Hamiltonian|77j . 
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Here, m, p, q and U {q, t) are the mass, momentum, position and potential of the particle, and 
and ujj are the mass, momentum, position and frequency variables of the jth bath 
oscillator mode. The quantities aj are coefficients that depend on the nature of the system-bath 
coupling. From eq ([T]), it is seen that the interaction part of the Hamiltonian is assumed to 
take the form Hi = —V{q)X, where V{q) is an arbitrary function of q and X = ajXj is 
the interaction coordinate. We introduced the counter term Ylj ^^q"^ /'^''^j^j maintain the 
translational symmetry of the Hamiltonian for U {q; 0) = 0. Here, we consider a reflection- 
symmetric potential system driven by a biharmonic force. [HI HSl [ISl [IZl IIHI UHl HSl HHl SZl SHI 
HH] This force prevents the system from reaching equilibrium and breaks the spatio-temporal 
symmetry. We write the potential as f/(g; t) = U{q) — qF(t), where U{q) is static, reflection- 
symmetric potential, and —qF{t) is the dynamic potential. Note that the time average of F{t) 
is zero. The total system obeys the von-Neumann equation (the quantum Liouville equation). 

After the bath degrees of freedom are traced out, the reduced density matrix elements of 
the system are obtained in the path integral form as 

p{<l,q'\t) = j dqodq'Q j Dq j Dq p{qo,qo)pcs{q, q' , t] qo,q'o) 

X exp{tSA[q,t]/h)F[q,q';t]exp{-iSA[q\t]/h), (2) 

where S'a[i?; t] is the action corresponding to the system's Hamiltonian, = /2m + U{q; t), 
p{qo, go) is the initial state of the system at time t^, and pcs{q, q', t; qo, q'^) is the initial corre- 
lation function between the system and the heat bath [78] . 

The effect of the bath is incorporated into the Feynman- Vernon influence functional, which 
is written 

g'; t] =exp|--^^*t^sV^(s) -^^ j\u'-^^{s - u)V°{u) + j\uC{s - u)V'' (u) j, 

(3) 



(1) 
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where Vit) = V{q{t)) - V{q'{t)) and V°{t) = V{q{t)) + V{q'{t)). Although the method we 
employ here can be used with any form of y(g).[59 | l60 | l6T] in this paper, we consider the linear- 
linear system bath coupling case, defined by V{q) = q. The interaction coordinate X = ajXj 
is regarded as a driving force through the interaction —qX. The canonical and symmetrized 
correlation functions, respectively, are then expressed as \E'(t) = f3{X; X{t))B and C(t) = 
X(0)})b, where (3 = l/ksT is the inverse temperature, X(t) is X in the Heisenberg 
representation, and (• • • )b represents the thermal average over the bath modes. |SH ES] Using 
the spectral density J{uj) = J2j cij^i^ — we can rewrite these functions as 

^(t) = 2 r cos(cjt), (4) 

Jo 

and 

POO 

C{t) = h I duo J {(jo) cos{(jjt) coth 
Jo 

The function C{t) is analogous to the classical correlation function of X{t) and corresponds to 
the correlation function of the bath-induced noise (fluctuations), whereas \E'(t) corresponds to 
dissipation. The function C(t) is related to \E'(t) through the quantum version of the fluctuation- 
dissipation theorem, C[uj\ = huJCoth{Phuj/2)/2'$[u)], which insures that the system evolves 
toward the thermal equilibrium state for finite temperatures, trB{exp[— /3Htot]}, in the case 
that there is no driving force |79] . 

As shown in Figure [T| the noise correlation C{t) becomes negative at low temperature, due 
to the contribution of the Matsubara frequency terms with u = 27r//3/i in the region of small t. 
This behavior is characteristic of quantum noise. The fact that the noise correlation takes neg- 
ative values introduces problems when the quantum Langevin equation is applied to quantum 
tunneling at low temperature. We note that the characteristic time scale over which we have 
C{t) < is determined by the temperature and is not influenced by the spectral distribution 
J{u)). Thus, the validity of the Markovian (or 5(t)-correlated) noise assumption is limited in the 
quantum case to the high temperature regime. Approaches employing the Markovian master 
equation and the Redfield equation, which are usually applied to systems possessing discretized 
energy states, ignore or simplify such tunneling contributions, often through use of the rotating 
wave approximation (RWA). The main problem created by the RWA is that a system treated 
under this approximation will not satisfy the fluctuation-dissipation theorem, and thus it may 
introduce significant error in the evolution of the system toward equilibrium. In the classical 
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Figure 1: The symmetric correlation C{t) defined by eq ([8]) is depicted as a function of tlie 
dimensionless time t for (a) tfie slow modulation case, 7 = 1, and (b) the fast modulation 
case, 7 = 10. The constants m(^7^/2 and h are set to unity. Note that 7 — )■ 00 corresponds to 
the Ohmic (Markovian) limit, as can be seen from eq (|6]). In each of the figures, the inverse 
temperatures are, from top to bottom, /3 = 0.2,0.5,1.0, and 5. The noise correlation C{t) 
becomes negative for (b), the fast modulation case, at low temperature (large (3) due to the 
contribution of Matsubara frequency terms. 



limit, with h tending to zero, C{t) is always positive. 

We assume that spectral density J {00) has an Ohmic form with a Lorentzian cutoff, i.e. 

m( 7^^ 



IT 7^ + 1^2' 



(6) 



where ( is the system-bath coupling strength, which represents the magnitude of damping, 
and 7 is the width of the spectral density of the bath mode. The canonical and symmetrized 
correlation functions then become 



-7|t| 



(7) 
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and 



cot 



,-"/\t\ 



(8) 



where V}^ = 2'Kk/f3h is the k-th. Matsubara frequency. |SU] Note that in the high temperature 
hmit, (3h'y <^ 1, the noise correlation function reduces to C{t) ~ m('je~"'^*^ / (3. This indicates 
that the heat bath oscillators interact with the system in the form of Gaussian-Markovian noise. 

The reduced hierarchy equations of motion (HEOM) can be obtained by considering the time 
derivative of the reduced density matrix elements given in eq g in the Wi gner representation. [56 
as a generalization of the Gaussian-Markovian quantum Fokker-Planck equation introduced by 
Tanimura and WolynesfSTf 158] . If we choose K so as to satisfy uk = 27!-K/{f3h) ^ Uc, where 
Uc is the characteristic frequency of the system, the factor e~'^'='*' in eq ^ can be replaced with 
the Dirac delta function, using the approximation i/fce~'^'='*'/2 ~ 5{t) (for k > K + 1), with 
negligible error. For the kernel equations ([T]) and Q, the number of hierarchy elements for 7 
is denoted by n, and the number of kth Matsubara frequencies is denoted by j'fc. 
the HEOM can be expressed asj^H E21 ES] 

K 

Cqm + S + n7 + ^ ikVk 



Then 



d 



k=l 
K 



k=l 



+ n7eo<":.']^(g,p;t) 



K 



(n) 



(9) 



fc=i 



where the quantum Liouvillian in the Wigner representation is expressed as [76] 



CQMW{q,p) = -^-W{q,p)- 
moq 



dp' 
2^ 



Uwiq,p-p';t)W{q,p'), 



(10) 



with Uw(yq,P',t) = 2 dr sm{pr / h)[U {q + r/2;t) — U{q — r/2;t)]. The other operators are 
defined as 



d_ 
dp' 



mh'j 
p H -— cot 



dp 



2m7^C d 
/3(z/| -72) dp' 

7 



(11) 
(12) 

(13) 



for > 0, and 

Note that the 0th element is identical to the Wigner function Wqq... oiQ^P'^t) = ^iQyP'jt), 
and the other elements are introduced in the numerical calculations in order to treat the non- 
perturbative, non-Markovian system-bath interaction. Although these elements do not have 
direct physical meaning, they allow us to take into account the quantum coherence and entan- 
glement between the system and the bath. [7H [75] The importance of the system-bath coherence 
was pointed out in the context of correlated initial conditions, [78] and it was shown with a non- 
linear response theory approach that the system-bath coherence plays an essential role in the 
case that the system is driven by a time-dependent external force. [56] 

The HEOM consist of an infinite number of equations, but it can be truncated at finite 
order with negligible error. [371 135] Essentially, the condition necessary for the error introduced 
by the truncation to be negligibly small is that the total number of hierarchy elements or the 
total number of Matsubara frequencies retained be sufficiently large. Explicitly, it can be shown 
that the condition N = n + Ylk=i3k ^ ^c/ min(7, z/i) is sufficient for this purpose. [81j For 
the Caldeira-Leggett Hamiltonian, the equations of motion are then truncated by using the 
" terminator" s, expressed in the Wigner representation as [6n 162] [63] 

i^it.A^^P-^t) = - {Cqm + e) (15) 

with n + Xlfc^i jk = N. The HEOM formalism can be used to treat a strong system-bath 
coupling non-perturbatively. It is ideal for studying quantum transport systems when employing 
the Wigner representation, because it allows the treatment of continuous systems utilizing open 
boundary conditions and periodic boundary conditions [76j. In addition, the formalism can 
accommodate the inclusion of an arbitrary time-dependent external field while still accounting 
for the system-bath coherence. Note that such coherence cannot be described if we assume that 
the state of the total system takes a factorized system-bath form. [56] Such system-bath coupling 
features are necessary to properly treat quantum ratchet systems. In the white noise (or 
Markovian) limit, 7 — > oo, which is taken after imposing the high temperature approximation, 
/3h'y ^ 1, the quantum Fokker-Planck equation can be derived in a similar form as the Kramers 
equation, [nU [52] which is identical to the quantum master equation without RWA.[82] Because 



1 cot 



K 

k=l 



Qp2 



(14) 
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we assume (3h'y -C 1 with 7 — )■ oo, this equation cannot be apphed to low-temperature systems, 
where quantum effects play a major role. 

The Wigner distribution function reduces to the classical one in the limit h ^ 0, and hence 
we can directly compare the quantum results with the classical results. The classical HEOM, 
which is derived from eq ^ by setting h = 0, is given by[571 158] 

{q,p;t)=- {CcL + n^) (g, p; t) 

+ n7c(p+'^i-]w^-'\q,p;t), (16) 



(3 dp J 



and the classical terminator is 
d 



l^(^) {q, p;t)=- {CcL + iV7) W^(^) {q, p; t) 



+ iV7c(p+^|)w^(^-^Hg,p;t). (17) 

The classical Liouvillian is defined by 

- Cc,W{q,p) = '^lw{q,p) + ^^lw{q,p). (18) 

The classical equation of motion is helpful, because knowing the classical limit allows us to 
identify the purely quantum mechanical effects [58^ l6T] . Note that the above equations of 
motion can also be derived from the classical Langevin equation, where the classical fluctua- 
tion dissipation theorem is hold.[5Z] Then, in the white noise (or Markovian) limit, the above 
equations reduce to the Kramers equation. 



In Sec. 4, we use the eqs (|9j)-(15) and eqs (16)-(18) to calculate quantum and classical 
ratchet currents, respectively. Because we are dealing with the distribution function for both 
quantum and classical cases, there is no need to sample kinetic trajectories, which is necessary in 
approaches employing the Langevin equation and in some kinetic approaches. This is convenient 
if we wish to calculate the current at low temperature, even in the classical case, because this 
sampling becomes inefficient due to the localization of the particle trajectories trapped in the 
potential. Another important aspect of the present methodology is that it allows elucidation of 
the dynamical behavior of the system through the time evolution of the distribution function in 
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the phase space. In addition to the dependence of the current on the temperature, the barrier 
height and the system-bath couphng, we analyze the mechanism of ratchet rectification using 
the Wigner and classical distribution functions. 



3 Potential, Force and Symmetry 

In the present study, we consider the one-dimensional symmetric potential (see Figure [2]) 

U{q) = UoCOs\Kq) (19) 

and the biharmonic force 

F{t) = Fi cos(fit) + F2 cos(21]t + 6), (20) 

where Uq is the barrier height, k is the wave number, and Fi, F2, VL and 9 denote the amplitudes, 
frequency, and the phase difference of the biharmonic forces, respectively. [HI UHl [161 [13 [IHl 
fT9| Hni Hni HH SHI Ho] Experimentally, such a situation has been reahzed for fluxons in long 
Josephson junctions[25] and for cold rubidium and cesium atoms in optical lattices. [261 EH 128] 
Before reporting the results of our calculations of the ratchet current, it is important to 
mention two symmetry conditions under which the current vanishes. [321 |33l [Ml [35l [36] The 
first one is the time-shift symmetry defined by the transformation {q,p,t) — )■ (— g, — p, t + T/2), 
where T is the period of the external force, and the second one is the time-reversal symmetry 
defined by the transformation {q,p,t) — > {q, —p, —t). When the equation of motion of the 
system is invariant under either of these two symmetry transformations, the current vanishes, 
because each of these transformations reverses the sign of the momentum. For the external 



force given by eq (20) and for the quantum and classical HEOM appearing in eqs (|9])-(15) 
and eqs (16)-(18), respectively, time-shift symmetry exists if Fi = or F2 = 0. Time-reversal 
symmetry is broken if 6* 7^ or tt in the Hamiltonian case, i.e. if C = 0, while, in the presence 
of dissipation, this symmetry is always broken, which implies that directed motion can exist 
for arbitrary 6. For the classical overdamped Langevin equation with white noise, the system 
possesses time-reversal symmetry if 6' = tt/2 or 37r/2. The overdamped (Smoluchowski) case 
can be studied within the HEOM formalism by choosing a large value for ( with 7 ^ in 
the classical case and by choosing a large value of ( with 7 ^ fi and (3h'j ^ 1 (i.e. the high 
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Figure 2: The Uq cos^(g) potential for the cases of (a) a low barrier, Uq = 1, (b) an intermediate 
barrier, Uo = 2, and (c) a high barrier, Uq = 4. In each case, eigenenergies and eigenfunctions 
with periodic boundary conditions are plotted as the black dashed lines and the red curves, 
respectively. 
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temperature limit) in the quantum case. In such cases, we expect the calculated current to 
vanish if 6 = 7r/2 or 37r/2. 



4 Numerical Results 

We numerically integrated the HEOM in the form of finite difference equations using the fourth- 
order Runge-Kutta method. We also imposed the periodical boundary condition W{p, q) = 
W{p, q+2n/K). The spatial derivative of the kinetic term in the Liouville operator, —{p/m)dW{p, q)/dq, 
was approximated using a third-order left-handed or right-handed difference scheme, depend- 
ing on the sign of the momentum, while other derivatives with respect to p were approximated 
using a third-order center difference scheme. [63] We employed the following units of the length, 
momentum, and frequency: qr = 1/ k, Pr = fin and Ur = hn? /m. Note that with these units, we 
have m = h = K, = 1. The dimensionless position and momentum then are given by q = q/qr 
and p = p/pr- The mesh sizes for the position and the momentum were chosen in the ranges 
0.0524 < Aq < 0.1047 and 0.04 < Ap < 0.10. The depth of the hierarchy and the number of 
Matsubara frequencies were chosen so as to satisfy G {5 — 10} and K E {1 — 2}. In this 
study, we fixed the inverse of the noise correlation time to 7 = l.Our and the amplitude and 
frequency of the biharmonic force to Fi = F2 = 0.20hu)r/qr and Q = l.Our, respectively. Three 
values were used for the height of potential barrier, Uq = huJr,2hujr, and Ahwr (see Figure |2]). 



Employing the biharmonic force given in eq (20), we integrated the classical and quantum 
HEOM until the distributions reached the steady state. Then, the classical and quantum 
directed currents (or the ratchet current) were evaluated from the distribution function as 



t+T 



J^lim^l dt'{p{t')), (21) 
where 

{p{t)) = j dqdppW{q,p;t). (22) 

We repeated the calculation for fixed physical conditions, while varying 6 in the range < 6 < 
27r. From these calculations, we find that as a function of 6, the current roughly takes the 
following form: [261 SH] 

J{9) = Jma.sm{e-eo), (23) 
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where Jmax and 6q are the maximum value of the current and the phase lag evaluated from the 
numerical simulations, respectively. Below we investigate the dependence of the current on the 
temperature, the barrier height and the bath coupling strength through the dependences of Jmax 
and 6q on these quantities. To gain further insight into the mechanism of ratchet rectification, 
we also present plots of the Wigner and classical distribution functions where helpful. 



4.1 Temperature Effects: the Classical Case 



Figure [3] displays (a) the maximum value, Jmax, and (b) the phase lag, 6*0, of the current as 
functions of the inverse temperature, (3, for three values of the barrier height, Uq, evaluated 



using the classical HEOM given in eqs (16)-(18). Here, we chose the weak coupling strength 



C = O.lOcUr- In all three cases depicted in Figure [sj^a), it is seen that the maximum current 
realizes a maximum value at some intermediate value of the inverse temperature. The value of 
P at which this maximum is realized is found to increase as Uq decreases. As shown in Figure 
[sj^b), similar profile is found for the phase lag 9q as a function of (3, although the peak positions 
are different. 

The behavior depicted in Figure |3] can be understood in terms of the thermal activation 
of a particle. To illustrate this, in Figure |4| we plot snapshots of the classical distribution 
function for (a) low, (b) intermediate, and (c) high temperature cases at five different values of 
the time for one cycle. In Figures |4]^a-i) and (b-i), it is seen that when the temperature takes a 
low or intermediate value, the distribution function is concentrated in the potential wells. The 
distributions are periodically accelerated by the biharmonic force and rotate clockwise around 
the bottom of each potential well. Figures |4]^a-ii) and (b-ii) reveal that when the distribution 
approaches the left side of the barrier, a small part of the distribution with positive momentum 
is transferred to the right potential well by crossing the barrier. This can be seen as the 
elongation of the distribution on the right side. In Figures |4|^a-iii) and (b-iii), it is seen that 
while the distribution begins to flow in the positive direction, the sign of the biharmonic force 
changes, and the main part of the distribution moves away from the barrier. The part of the 
distribution that crosses the barrier flows into the right potential and rotates clockwise in a 
spiral form while losing energy, due to dissipation. Figures |4]^a-iv) and (b-iv) show the flow 
stops shortly after the biharmonic force change direction, with a small time delay. In Figures 
|4|^a-v) and (b-v), it is seen that when the distribution approaches the left-hand side of the 
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Classical 



Quantum 




Figure 3: The maximum current in (a) the classical and (a') the quantum case and the cor- 
responding phase lag in (b) the classical and (b') the quantum functions of the in- 
verse temperature, /3, for three values of the barrier height Uq. The coupling strength here is 
C = O.lOw^. 
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Figure 4: Snapshots of the classical distribution for 6 = O.Tvr in the case of intermediate barrier 
height, Uq = 2, for three values of the inverse temperature: (a) low temperature, /3 = 3; (b) 
intermediate temperature, /3 = 2; (c) high temperature, 13 = 1. The snapshots correspond to 
the following times: t = (i) 0, (ii) OAn, (iii) O.Svr, (iv) 1.27r, (v) l.Gvr (l/Q). 
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barrier, a part of the distribution with negative momentum also flows in the opposite direction, 
but the effect of the acceleration caused by the biharmonic force is smaller in this case than 
in the positive case, due to the asymmetric influence of the biharmonic force on the system 
dynamics with the phase lag Oq. 

In the intermediate temperature case depicted in Figure |4]^b), the distribution is broadened 
compared with Figure |4]^a) due to thermal effects. The current is greater than in the low 
temperature case, because the weight of the distribution at energies above the barrier under the 
biharmonic force is greater. The response of the distribution to the external force is also delayed 
due to the broadening of the distribution in the p direction. Because the phase lag, 60, changes 
in accordance with the timing of the excitation created by Fi cos{Qt)q and F2 cos{2Qt + 6 — 6Q)q 
through the gradients of the potential, this delay may be the reason that ^0 increases as (3 
decreases up to the maximum point. 

In the high temperature case considered in Figure |4] (c), the distribution is populated even 
at the top of barrier. This dispersion of the distribution indicates that the effect of the ratchet 
potential is suppressed at high temperatures. This causes the spatio-temporal asymmetry, 
which is necessary for ratchet current, to be weak, and as a result, the current decreases as 
the temperature increases. Since the activation energy increases as the barrier increases, the 
maxima of the peaks shift toward higher temperature as Uq increases. In the high temperature 
case, because the distribution is dispersed over the potential, the distribution can flow to the 
neighboring potentials both in the positive and negative directions with only a small delay. This 
can be understood from Figure |4]^c-ii), where it is seen that the current flows quickly following 
the movement of the distribution. Thus, 60 decreases as the temperature increases, although 
the peak positions are different from those of Jmax- This is due to the fact that there is only 
an indirect relation between the current and the phase. 



4.2 Temperature Effects: the Quantum Case 

Next, we consider the quantum mechanical case under the same physical conditions as in the 
classical case considered above. Figure [s] displays (a') the maximum value, Jmax, and (b') the 
phase lags, of the current as functions of the inverse temperature, /3, for three values of the 
barrier height, Uq, calculated using the quantum mechanical HEOM given in eqs (|9])-(15). It is 



seen that the profiles of Jmax and 60 differ significantly from those in the classical case. In the 
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intermediate and high barrier cases, The figure also shows that the quantum mechanical Jmax is 
significantly larger than the classical one for the cases of an intermediate and high barrier in the 
low temperature regime (/3 > 2), while they exhibit similar behavior in the high temperature 
regime (/3 < 1). 

To see a role of quantum effects, we consider the crossover temperature at which the classical 
thermal-activated regime crosses over to the quantum tunneling regime with regard to chemical 
reaction rates |H3] for the bath spectral density, J (to), given in eq ^ with the potential heights 
Uq = 1, 2, and 4. The evaluated crossover temperatures are at /3c=4.5, 3.2, and 2.2, respectively. 
These values of Pc are much smaller than those at which we observed the difference between the 
quantum and classical cases. This suggests that our analysis based on the reaction rate may 
not apply directory, because here we considered Jmax under the external driving force while 
altering the phase. 

A distinctive feature of the quantum mechanical results is found in the low barrier case (the 
red curve) in Figure [sj^a'). One may expect that the ratchet current in the quantum mechanical 
case is larger than that in the classical case, especially at low temperatures due to tunneling. 
However, our results in the low barrier case presented in Figure [3|^a') reveal that in fact this is 
not the case. 

To elucidate the reason for the behavior depicted in Figure |3]^a') for the case Uq = 1, we 
consider the Wigner distribution in the low and intermediate barrier cases for three values of the 
temperature in Figures |5] and [6| respectively. In both cases, the Wigner distribution function 
is dispersed over the potential, and it is negative in some regions (the regions inside the white 
and black loops) at low temperature, (/3 > 2). The Wigner function for the quantum system is 
not positive definite, and in cases that the tunneling process is important, the region in which 
it takes negative values is larger. [Mj This indicates that delocalization of the distribution in 
the low barrier case depicted in Fig [s] (a) arises from tunneling. This can be understood from 
Figure [2]^a), where the ground state in the low potential case is seen to have a large population 
under the barrier. The tunneling contribution induced by biharmonic perturbation may be 
suppressed at high temperature, because the temperature enters the theory through f3h, and 
hence the classical limit, h 0, is equivalent to the high temperature limit, /3 — 0. In this 
case, however, thermal activation causing fiow over the potential barrier is significant. For this 
reason, the distribution is dispersed for all temperatures. As shown in the high temperature 
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Figure 5: Snapshots of the quantum distribution for 9 = OAtc in the case of low barrier height, 
Uq = 1, calculated using the HEOM for three values of the inverse temperature: (a) low 
temperature, (3 = 3; (b) intermediate temperature, (3 = 2; (c) high temperature, [3 = 1. These 
snapshots correspond to the following times: t = (i) 0, (ii) OAtt, (iii) O.Svr, (iv) 1.2n, (v) l.Gvr 

(i/n). 
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Figure 6: Snapshots of the quantum distribution for 9 = O.Gvr in the case of intermediate barrier 
height, Uq = 2, for the following values of (3: (a) low temperature, (3 = 3; (b) intermediate 
temperature, /3 = 2; (c) high temperature, /3 = 1. These snapshots correspond to the following 
times: t = (i) 0, (ii) OAtt, (iii) O.Stt, (iv) 1.27r, (v) I.Qtt 
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classical case depicted in Figure |4[b), the effect of the ratchet potential is weak when the 
distribution function is dispersed. This is the reason that the ratchet current is smaller in the 
quantum case than the classical case if the barrier is low. 

In the cases of intermediate and high barriers, depicted in Figure [3} tunneling also plays 
a significant role. As shown in Figure [6} the Wigner distribution is not dispersed when the 
barrier is high. As a result, the ratchet mechanism is effective with the help of tunneling, and 
thus the quantum mechanical current is larger than the classical current. However, when the 
temperature becomes very high, the quantum coherence is lost due to the thermal noise, and 
as a result, the quantum value approaches the classical one. 

As shown in Figure |3|^b'), the phase lag has a critical point in the high temperature region 
(/3 < 1.5) as in the classical case depicted in Figure [3|b). However, in the low temperature 
region(/3 > 2.5), the phase is seen to increase monotonically in the quantum mechanical case, 
while it decreases monotonically in the classical case. We carefully investigated the dynamics 
of the Wigner distribution in the low temperature region by studying snapshots of it and the 
corresponding current using fine time slices. We thus found that the transfer of the distribution 
through tunneling takes a longer time at low temperature. This is because thermal transitions 
between the ground and excited states, which is the cause of tunneling transition induced by 
biharmonic perturbation, may be suppressed at low temperature. Because the net current flows 
slowly in the low temperature case, the phase lag that controls the timing of the biharmonic 
force increases when the temperature decreases. For higher barriers, the tunneling plays a 
minor role, and thus the increase of the phase is smaller than in the low barrier case. 

Note that, as discussed in Ref. [85], quantum diffusion may be suppressed by quantum 
reflection. In the cases we studied, however, the curvature at the top of the barrier and in 
the potential well are equal, and the force that creates the difference between the curvatures is 
small. [86] Thus, in the present case, quantum reflection is not the cause of the suppression. 

4.3 Varying the Coupling Strength 

We next investigate the change of the ratchet current as the system-bath coupling strength ( 
is varied by taking advantage of the fact that the HEOM constitute a non-perturbative theory. 
In Figure [7| we plot the classical (red curve) and quantum (green curve) ratchet current for 
(a) the low barrier case, Uq = 1, and (b) the intermediate barrier case, Uo = 2, with (3 = 1. 
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Figure 7: The maximum ratchet current in the classical (red curve) and quantum mechanical 
(green curve) cases for (a) a low barrier, Uq = 1, and (b) an intermediate barrier , Uq = 2, with 
(3 = 1. The insets depict the phase lag 9q as a function of (. 
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(a)C = 0.1 (b)C = 0.3 (c)C = 0.5 (d)C = 0.7 
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Figure 8: Snapshots of the classical (upper panel) and quantum mechanical (lower panel) 
distribution functions at the fixed time t = 2tt (1/f^) for four values of the coupling strength: 
C = (a) 0.1, (b) 0.3, (c) 0.5, (d) 0.7. The upper panels, (i) and (ii), describe the classical cases, 
whereas the lower panels, (i') and (ii'), describe the quantum cases for the case of a low barrier, 
Uo = I, and an intermediate barrier, Uq = 2, respectively. 



22 



In both cases, the classical and quantum currents seem to approach the same value for large 
although the quantum values reach the large ( limit more rapidly. To compare the two 
cases, we plot the quantum and classical distribution functions in Figure |8] It is seen that in 
the low barrier case, the quantum distribution is dispersed, as in the case depicted in Figure 
[Sj Moreover, in the quantum mechanical case, the distribution in the p direction becomes a 
symmetric Gaussian for ( > 0.3, while there is some deviation from Gaussian in the classical 
case, due to the external perturbation. This Gaussian feature arises from quantum tunneling, 
as can be observed more clearly in the low barrier case. It is interesting that the situation in 
which there exists a symmetric Gaussian distribution in momentum space is similar to that 
described by the classical Smoluchowski equation, which can be derived in the overdamped 
limit from the Kramers equation. In the overdamped limit, the distribution of the momentum 
direction is Gaussian. In the presently considered case, we find that such a situation arises due 
to the effect of tunneling even if the system-bath coupling is weak. Therefore, it is natural to 
call this situation a tunneling-induced Smoluchowski limit. 

The phase approaches n/2 as C increases, as depicted in the inset of Figure [7](ii). This is 
similar to the situation for the classical Smoluchowski equation. For Uq = 2, the quantum 
and classical cases exhibit similar ( dependence. This is due to the fact that the dispersion 
of the distribution through tunneling does not occur in the higher barrier depicted in 

Figure [6j The quantum values approach the classical values for large (, because in that case, 
the quantum coherence is destroyed by the system-bath coupling. 

5 Concluding Remarks 

The classical and quantum mechanical ratchet currents in a system driven by a biharmonic force 
were calculated in a rigorous manner for the first time using the reduced hierarchy equations of 
motion (HEOM) in the Wigner representation over a wide range of values of the temperature, 
the system-bath coupling strength, and the potential height. The roles of fluctuations, dissipa- 
tion, and the biharmonic force on the ratchet current were investigated in classical and quantum 
mechanical cases by studying the dependence of the ratchet current on these parameters. 

In the classical case, we found that over the temperature range we studied, the current 
realizes a maximum at an intermediate value, and falls off for both high and low temperatures. 
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This decrease on the high temperature side is a result of thermal activation, which induces 
dispersion of the distribution over the barrier in the high temperature region. It the quantum 
mechanical case, that tunneling enhances the current in the high barrier case, while it suppresses 
the current in the low barrier case. This is because the effect of the ratchet potential is weak 
in the case of a low barrier, due to the dispersion of the distribution that arises from tunneling. 
Moreover, in the quantum mechanical case, there is a classical Smoluchowski-like regime induced 
by tunneling effects. We also found that the phase lag of the biharmonic force is related to the 
response time of the system to the external perturbation. 

As we demonstrated, the dynamical behavior of the system is clearly and readily elucidated 
by the time evolution of the Wigner distribution function. Although we must introduce hier- 
archies of Wigner functions described by a discretized mesh, a modern personal computer is 
sufficiently powerful to solve the equations of motion in a reasonable amount of time. Numer- 
ical techniques, such as the optimization of the hierarchy. [871 188| \89\ [90] the utilization of a 
graphic processing unit (GPU),[91J and memory allocation for massive parallel computing P2] 
have been developed for the HEOM approach in order to accelerate numerical calculations. As 
an approximated method, it may be possible to adapt a multi-configuration time-dependent 
Hartree (MCTDH) approach, which is shown to be an efficient method to solve Schrodinger 
equations of motion. [931 IMl 195] With these developments, the extension of the present study to 
the variety of multi-dimensional potential systems with non-Drude-type spectral distribution 
functions may be possible. 
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